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a&stract 


Neu  labeling  technique,  are  provided  for  accelerating  the  basis 
exchange  step  o,  specialised  linear  programing  methods  for  netuorh 
problems.  These  technics  substantial!,  rednce  the  amooht  of  computation 

involved  in  updating  operations. 


1.0  INTRODUCTION 

In  solving  mini™  cost  flow  network  problems  by  specialized  linear 
programming  methods,  an  important  question  is:  how  can  one  update  the 

ehl.  question  is  provided  by  special  list  structure  techniques  such  as  the 

_ ATT  _.^hocl  r a 1 which  have  contributed 
API  method  [11  and  the  more  recent  ATI  method 

. .no  rhp  efficiency  of  network  algorithms  (see,  e.g., 

dramatically  to  improving  the  efficiency 

[2,  3,  6,  7,  8!).  This  paper  addresses  the  issue  o£  which  supplemental 
techniques  can  be  used  to  enable  these  list  structures  (and  particularly 

the  ATI  method)  to  be  implemented  with  greater  efficiency. 

As  shown  in  [1],  the  major  updating  calculations  of  a basis  exchange 

seep  can  be  restricted  to  just  one  of  the  two  subtrees  created  by  dropping 

i 4c  rn  identify  the  smslleir 

the  outgoing  arc.  Consequently,  a natural  goal  is 

nf  . i q of  c(x)  that  names  the  number  of  nodes 
of  these  two  subtrees  by  means  of  a list  t(x)  tna 

. v „ A r lever  and  rather  intricate  procedure 
in  the  subtree  "headed  by  node  x . A clever  ana 

for  doing  this  -as  proposed  by  Srinlvasan  and  Thompson  19).  However,  unfor- 
tunately, this  procedure  requires  sorting  the  nodes  of  the  subtree  by  their 
distances  from  the  toot,  and  then  further  entails  a full  subtree  update  of 
both  the  distance  values  and  the  t(x>  values  at  each  basis  exchange  step. 

Because  of  the  substantial  amount  of  work  required  to  update  the  t(x)  list, 

its  use  has  never  been  found  practicable. 

The  purpose  of  this  paper  is  to  propose  a new  type  of  relabeling  scheme 

that  succeeds  in  updating  t(x>  without  sorting.  In  fact,  this  scheme  requires 
even  less  work  than  to  update  the  distance  values  of  ,9).  The  relabeling  is 
based  on  "absorbing"  C(x>  into  the  updating  calculations  of  the  ATI  method. 
Moreover,  these  calculations  are  carried  out  simultaneously  with  the  procedures 
(4)  for  updating  other  changes  introduced  by  the  basis  exchange  step. 
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To  achieve  the  integration  of  the  ATI  calculations  and  the  update  of 
t(x),  we  introduce  an  index  function  f(x)  that  names  the  last  node  in  the 
subtree  rooted  at  x.  We  show  that  f(x)  additionally  makes  it  possible  to 
streamline  the  ATI  calculations.  Finally,  as  a bonus,  we  show  that  t(x) 
can  accommodate  all  of  the  relevant  functions  filled  by  the  distance  values, 
and  hence  can  replace  these  values.  The  net  gains  of  all  these  advantages 
produce  a substantially  improved  procedure  for  implementing  the  basis 


exchange  operations. 
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2.0  NOTATION  AND  DEFINITIONS 

We  shall  follow  the  notational  conventions  of  [1],  identifying  the  components 

of  the  basis  exchange  step  as  follows: 

(p , q)  = the  arc  leaving  the  basis,  where  p is  currently  the 
predecessor  of  q (via  the  'predecessor  indexing"  used  with  the  ATI  method) . 

(u,v)  = the  arc  entering  the  basis,  where  u is  the  node  whose 
"backward  path"  (consisting  of  all  its  ancestors  under  the  predecessor 
indexing)  contains  arc  (p,q). 

T = the  basis  tree 

T(x)  = the  subtree  of  T that  is  rooted  at  node  x (hence  the  sub- 
tree that  includes  x and  all  its  descendants  under  the  predecessor  indexing) 
S(x)  = the  "thread  successor"  of  x as  defined  in  the  ATI  method 
[A],  (In  particular  T(x)  = {s°(x),  s^(x),  ....  sr(x)}  for  some  r > 0, 

where  s°(x)  = x and  s^+^(x)  = s(s^(x)).) 

f(x)  = the  "last  node"  ("sr(x)")  in  T(x) 
t(x)  = the  number  of  nodes  in  T(x) 

The  basis  exchange  step  may  be  visualized  as  consising  of  two  components: 

(1)  Dropping  arc  (p,q)  to  create  two  Independent  subtrees:  T(q) 

and  T - T(q)  (where  the  latter  is  the  subtree  of  T that  excludes  T(q)  and 

all  its  nodes,  and  hence  which  excludes  the  "connecting  arc"  (p,q)); 

(2)  Adding  arc  (u,v)  to  create  a single  new  basis  tree. 

The  subtrees  T(q)  and  T - T(q)  can  be  viewed  as  any  two  node-disjoint  trees 

which  are  to  be  joined  by  an  arc  to  create  a new  basis  tree.  We  shall  therefore 

first  develop  updating  operations  to  make  T(q)  and  T - T(q)  into  label  "indepen- 
dent" trees  in  preparation  for  selecting  which  is  to  be  the  new  "upper  tree" 
(called  T^)  rooted  at  x ^ and  which  is  to  be  the  new  "lower"  trees  (ca'.lec  7.) 
rooted  at  We  may  then  assume  that  the  root  of  T becomes  the  root  of  the 

new  basis  tree.  Additionally,  we  let  (y^,  y be  the  arc  that  joins  T^,  and  T2 
where  y^  is  a node  of  T^  and  y2  is  a node  of  T2.  NSxt  we  develop 
operations  to  re-root  T2  at  y2  in  preparation  for  attaching  T2  to  T-^  via  arc 
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( yp  y 2)  to  create  the  new  basis  tree.  (There  is  no  requirement  that  be  dis- 
tinct from  x^  or  that  yj  be  distinct  from  X2«) 


3.0  UPDATING  OPERATIONS 

Using  the  preceding  definitions,  we  show  how  to  find  updated  values  s*(x), 
t*(x) , f*(x)  for  s(x),  t(x)  and  f (x)  as  follows. 

I.  Update  s(x),  t(x),  f(x)  to  make  the  subtrees  T(q)  and  T - T(q)  independent 


1.1.  Update  for  T - T(q): 

For  s*(x):  Identify  the  node  y in  T - T(q)  such  that  s(y)  * q. 

Then  set  s*(y)  = s(f(q)).  No  other  s(x)  values  are  changed.  (Note:  while 
this  step  is  obviously  facilitated  by  the  direct  accessing  f(q),  the  identi- 
fication of  y is  further  speeded  by  utilizing  the  f(x)  functions  as  follows: 

First,  let  y1  = p.  Second,  if  s(y^)  = q,  then  y = y^  and  the  process  stops. 

•1  1 

Otherwise,  let  y = f(s(y  ))  and  repeat  the  second  step.) 

For  t*(x):  Set  t*(x)  = t(x)  - t(q)  for  those  nodes  x on  che  "back-  • 

ward  path"  from  p to  the  root  of  T.  Due  to  cancellation  effects  of  subse- 
quent calculations,  this  step  can  be  restricted  to  the  partial  backward  path 
from  p to  the  "intersection"  node  z that  is  the  unique  node  of  the  basis 
loop  that  lies  on  the  backward  paths  from  both  u and  v,  excluding  node  z 
itself  from  consideration.  (Thus,  possibly  there  may  be  no  updating  in  this 

step . ) 

For  f “ (x)  : Set  f''(x)  = y for  those  nodes  x on  the 

predecessor  from  p to  the  root  such  that  f(x)  = f(q).  (If 
f(p)  i f(q)  no  updating  is  done.) 

1.2.  Update  for  T(q) : 

For  s*(x) : Set  s*(f(q))  = q. 

No  other  updating  of  any  s(x),  t(x)  or  f(x)  is  required  for  T(q} 


II.  Decide  which  of  T(q)  and  T - T(q)  is  to  be  T,  and  which  is  to  be  T„. 
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„ t(q)  exceeds  naif  the  number  of  nodes  In  the  network  let  TCq)  be  V and 
let  T - T(q)  be  T,  (hence  x£  - , and  x,  - the  root  of  T) . Otherwise,  let  T - T(q) 

be  T , and  let  T(q)  be  T2* 

The  principal  reason  for  re-rooting  Is  to  minimise  the  computational  effort 
when  the  updating  of  the  dual  variables  is  integrated  with  the  other  updating 
operations.  However,  it  is  possible  to  separate  these  updating  operations,  carryin, 

out  the  node  potential  update  on  the  smaller  subtree,  and  always  earning  out  the 

aihfrpe  T(fl)  thereby  leaving  the  root  unchanged, 
other  updating  operations  on  subtree  (q;, 

f niu,  directly  from  the  subsequent  development. 

This  latter  case  follows  directly 

, T reversing  the  predecessor  orientation  of  the  path 
III.  Make  y2  the  new  root  of  T2,  reversing 

from  yo  t0  x2# 

, , r /v\  -i q parried  out  exactly  us 

For  s*(x):  This  portion  of  the  updating  o 

specified  in  [4]. 

/ h /r,  j f Tfni  = t , and  the  restricted  update 

For  t*(x):  Let  t*(y,)  - t(x2).  (But  if  T(q)  Tj. 

of  t (x)~wascarr  led  out  in  step  1.1  - which  Is  computationally  preferable  - set 

t*(y  ) . t(x  ) - t(q>.)  Then  for  each  x on  the  path  from  x.,  to  y2>  exclu  ing  x y2 

a t*(x)  - t*(y  ) - t(x),  where  5 is  the  successor  of  x via  the  predecessor  index!  .g 

S C C 2 . fTt  is  important  to  note  here  that 

before  reversing  the  orientation  of  the  path,  (It  imp  ltlnttte 

t(x)  and  t*(x)  must  be  kept  distinct  from  each  other;  I.e., 

to  replace  «(S)  by  t*U)  before  computing  ftt I = t*  y ) - « ^ ^ 

For  f*(x) : The  node  of  T2  that  is  examined  last  m P 

, u T via  the  ATI  method  is  designated  to  be  f*(x2).  n 
the  node  potentials  2 £(.  ^ „here  - ls  the  successor  of  X,  before 

particular,  if  x2  y2,  or  2 2 (i.e.,  no  change  occurs 

reversing  the  orientation  of  the  path,  then  £ *2  ’ (This 

, )co  f*lx  ) = v for  the  node  y such  that  s(y)  2 
in  f(xJ.)  Otherwise,  f U2'  y , conciusion 

2C  . . , i but  it  „ill  be  identified  automatically  at  the  concx 

may  be  found  as  in  1.1,  but  wll  all  x on  the 

of  updating  node  potentials.)  Thereupon,  set  f <x>  2> 

path  from  *2  to  ,r  (No  other  changes  are  made.) 

/ \ f he  new  basis  tras 

IV.  At-t-ach  T:  to  T]  

(,,h*re>  T2  is  now  rooted  at  as  a result  of  steB_m).. 

For  s*(x).:  Set  s*(f*(y2»  = S<V  and  s*(yl}  = Y2* 
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Fur  t*(x) : Set  t*(x)  = t (x) + t* (y 2)  for  all  x on  the  path  from 

y ^ to  Xj.  (But  if  T(q)  = T?,  and  the  restricted  updaLe  of  t(x)  was  applied 
in  step  £.1,  then  the  current  step  should  be  restricted  to  those  x on  the 
path  from  y^  to  z,  excluding  z itself,  for  the  "intersection"  node  z as 

identified  in  1.1.)  __  

For  l*(x):  If  fty,)  t V sot  f*(x)  = f<y2>  for  those  nodes 

x on  the  predecessor  from  y-  to  x,  such  that  f(x)  - fty^)-  OJ 5 

changes  are  made  if  f(y  ) r y ^*) 

The  proposed  procedure  for  updating  t(x)  clearly  requires  less  effort  than 
updating  the  distance  function  of  [9],  which  minimally  involves  an  addition  for 
every  node  of  the  subtree  T(q),  and  in  the  case  of  T(q)  = T^,  requires  an  addition 
for  every  node  of  T.  The  fact  that  t(x)  can  replace  the  distance  function  is  a 
direct  consequence  of  the  observation  that  if  t(x)  < t(x^),  then  x cannot  be  a 
descendant  of  x*.  Thus,  t(x)  can  be  used  in  essentially  the  same  manner  as  the 
distance  function  to  facilitate  operations  such  as  involved  in  identifying  the 
loop  created  by  adding  arc  (u,v)  to  the  basis  tree. 

4.0  INITIALIZATION 

It  is  left  to  characterize  the  procedure  for  establishing  the  initial  values 
of  t(x)  and  f(x).  This  occurs  simultaneously  with  the  initial  determination  of 
the  s(x)  values  as  follows. 

Consider  the  step  in  which  s^+^  (x0)  = s(s^(x0))  is  identified  (K  > 0). 

If  b^(xo)  is  the  predecessor  of  s^+^(x0)  (via  the  predecessor  indexing),  do 

j 

nothing.  Otherwise,  for  all  nodes  s1(xe)  on  the  backward  path  from  s (x0)  tc 

k+1  k+1 

the  predecessor  of  s (x0),  excluding  the  predecessor  of  s (x0)  itself,  set 

t(8^(x0))  = k+1  - i,  and  set  f(s*(x0))  a s^(x0). 

When  the  last  node  sn  ^(x0)  of  the  network  is  determined  (where  n = the  total 

number  of  nodes  in  the  network),  set  t(s^(x0))  = n - i and  set  f(s^(x0))  = 

n_  ^ 2 n~  1 

s (x0)  for  all  s (x0)  on  the  backward  path  from  s (x0)  to  xe. 
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To  easily  keep  track  of  Che  index  1 for  each  node  s1^.)  that  Is  to  be 
considered  on  a given  step,  it  is  convenient  to  keep  a list  that  consists  pre- 
cisely of  the  indexes  i of  the  nodes  .*(..)  to  x„ . Specifically,  to  begin  the 
list  contains  the  single  index  0 (for  s‘(x.».  When  s*+1<x.)  is  created,  the 
number  k+l  is  added  to  the  end  of  the  list.  When  a backward  path  from  s <*.) 
is  traced,  consisting  of  r nodes  (say)  s‘(x.)  whose  values  t(s\x„»  and  f(s.(x„» 
are  to  be  set,  the  indexes  of  these  r nodes  will  be  exactly  the  corresponding 
last  r numbers  on  the  list.  By  removing  these  numbers  from  the  list  just  before 
adding  the  number  k+l,  the  desired  structure  of  the  list  is  maintained. 

5.0  CONCLUDING  REMARKS 

The  computational  advances  in  network  algorithms  afforded  by  the  techniques  to 
[1,  4,  5,  9], and  the  evident  improvements  in  these  techniques  provided  by 
forgoing  procedures,  argues  strongly  in  favor  of  the  implementation  value  of  these 
procedures.  Moreover,  we  conjecture  that  coupling  the  ATI  method  with  the  use 
of  t(x)  and  f(x)  in  the  manner  described,  not  only  dominates  all  previous  updating 
procedures,  but  closely  approaches  the  "asymptotic  limit"  of  effective  trade-off 
between  computational  efficiency  and  memory  requirement  for  a broad  class  of  net- 
works. Results  of  a computational  study  of  these  new  procedures,  imbedded  in 
various  alternative  Implementations  of  specialised  simplex  codes,  will  be  reported 


in  a subsequent  paper. 
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